Astronomy & Astrophysics manuscript no. vilhu 


© ESO 2009 


April 25, 2009 





o 
o 

< 

(N 



Oh 
6 

o3 



> 

On 

m 

O 

On 
O 



13 



Orbital modulation of X-ray emission lines in Cygnus X-3. 

0. Vilhu 1 , P. Hakala 2 , D.C. Hannikainen 1 - 4 , M. McCollough 3 , and K. Koljonen 3 ' 1 

1 Observatory,Tahtitominmaki (PO Box 14), FI-00014 University of Helsinki, Finland 
e-mail: osmi . vilhu@astro . helsinki . f i 

2 Tuorla Observatory, Vaisalantie 20, FI-21500 Piikkio, Finland 
e-mail: pahakala@astro.utu.fi 

3 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 
e-mail: mccolml@head-cf a . harvard . edu 

4 Metsahovi Radio Observatory/TKK, Metsahovintie 1 14, FI-02540 Kylmala, Finland 
e-mail: diana@kurp . hut . f i 

Received 05/1 1/2008; accepted 19/04/2009 

ABSTRACT 

Aims. We address the problem where the X-ray emission lines are formed and investigate orbital dynamics using Chandra HETG 
observations, photoionizing calculations and numerical wind-particle simulations. The aims were to set constraints on the masses of 
the components of this close binary system consisting of a Wolf-Rayet (WR) star and a compact component and to investigate the 
nature of the latter (neutron star or black hole). The goal was also to investigate P Cygni signatures in line profiles. 
Methods. The observed Si xiv (6. 185 A) and S xvi (4.733 A) line profiles at four orbital phases were fitted with P Cygni-type profiles 
consisting of an emission and a blue-shifted absorption component. Numerical models were constructed using photoionizing calcula- 
tions and particle simulations. In the models, the emission originates in the photoionized wind of the WR companion illuminated by 
a hybrid source: the X-ray radiation of the compact star and the photospheric EUV-radiation from the WR star. 
Results. Spectral lines with moderate excitation (such as Si xiv and S xvi) arise in the photoionized wind. The emission component 
exhibits maximum blue-shift at phase 0.5 (when the compact star is in front), while the velocity of the absorption component is con- 
stant (around -900 km/s). Both components, like the continuum flux, have intensity maxima around phase 0.5. The simulated Fexxvi 
Lya line (1.78 A, H-like) from the wind is weak compared to the observed one. We suggest that it originates in the vicinity of the 
compact star, with a maximum blue shift at phase 0.25 (compact star approaching). By combining the mass function derived with that 
from the infrared He i absorption (arising from the WR companion), we constrain the masses and the inclination of the system. 
Conclusions. The Si xiv and S xvi lines and their radial velocity curves can be understood in the framework of a photoionized wind 
involving a hybrid ionizer. Constraints on the compact star mass and orbital inclination (i) are given using the mass functions derived 
from the Fe xxvi line and He i 2.06 yum absorption. Both a neutron star at large inclination (;' > 60 degrees) and a black hole at small 
inclination are possible solutions. The radial velocity amplitude of the He n 2.09 /im emission (formed in the X-ray shadow behind 
the WR star) suggests i = 30 degrees, implying a possible compact star mass between 2.8 - 8.0 M©. For i = 60 degrees the same 
range is 1.0-3.2% 

Key words. Stars:individual:Cyg X-3 - Stars:binaries:spectroscopic - Stars:winds - Accretion - Stars:neutron - Black hole physics 

1. Introduction spectrum (Paerels [2000l ). The spectrum shows radiative recom- 

^ „ .„ „ s . ,. , bination continuums (RRC) that is a clear indication of a pho- 

^S^iis X-3C4U 2030+40, V1521 Cyg) is a high-mass X-ray toionized lasma (an altern ative to a collisionally ionized gas), 

binary (HMXB) located atadistance of 9kpc with a close binary Uedahl et aJ {2Qm md McColl h et aL ^000) report P 

orbit (P = 4.8 hour; Liu|2002l Hanson 2000) The compact star Q { fileg ^ some x _ emission lines . Stark and Saia 

is either a neutron star or a black hole and the companion is a ^ found D j modulation of X -ray lines in Chandra 

WN5-7 type Wolf-Rayet (WR) star (van Kerkwijk| 1992|199^ , A ^ when c x _ 3 was M h and 

with strong evidence pointing to an early-type WN star (Fender used Ms tQ constrain the compact star mass to below 3 6 solar 



masses. 



et al. 1999). The system appears to be engulfed in an extended 

envelope produced by the wind of the WR star. This notion is 

Z j , c At • j 11 i_ j In the present paper we analyse one long Chandra HETG 

supported by the tact that massive winds are generally observed . . , \ . .. , . _ „ ^ , 

in Wolf-Rayet stars (Langer£989, Crowther 2007). observation taken during a high state of Cyg X-3. We fit the 

Specifically with regard to Cyg X-3, large and complex ab- Si XIV 6.185 A Lycr-line at four orbital phases with a P Cygni- 

sorption is found in br oadban d X-ray spectral studies (e.g. Vilhu ^ P rofile < an eimsslon component with a blue-shifted absorp- 

120031 Hialmarsdotter l2008a1 Szostek and Zdziarski 120081 , the tlon component). Results for the S xvi Lya-line at 4.73 A, the 

orbital modulation of X-ray light curves can be attributed to Fe Lya line at 1.78 A, and the Hea line at 1.859 A are also pre- 

electron scattering absorption in a n asym metric wi nd (wi th re- sented. Using a two-component hybrid ionizing flux arising from 

spect to the X-ray emitter, Pringle \WT4\ Willingale [T983T l, and *e compact object a nd the WR star and the XSTAR photoion- 

the source exhibits a very rich Chandra HETG emission line izing code (KallmanJgXll for a review of photoionized winds 

see also Liedahl 120011 ). we estimate line emissivities in the ra- 

Send offprint requests to: O. Vilhu diation field. Radial velocities and lightcurves were computed 
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from a wind model based on particle simulations in the binary 
potential and compared with those we observed. We conclude 
by discussing constraints on the orbital inclination, component 
masses, and wind velocity. 

2. Chandra observations 

Chandra High Energy Transmission Grating (HETG) observa- 
tions (OBSID: 7268, PI: M. McCollough) were scheduled on 
2006 January 25-26 as a part of a large international campaign 
to study Cyg X-3 during an active state. The observations started 
on MJD 53760.59458 (when the source was at phase 0.053) and 
went through to MJD 53761.42972 (phase 4.240). During these 
observations Cyg X-3 was in a high soft state (quenched radio 
state) with an average RXTE/ ASM count rate around 30 cps (cor- 
responding to 400 mCrab; typical hard state count rates are less 
than 10 cps). 

The HETG consists of two gratings, the Medium Energy 
Grating (MEG) and the High Energy Grating (HEG). The ob- 
servations were done in CC mode (continuous clocking) with 
a window filter applied to the zeroth order. The CC mode was 
used to avoid pileup in the dispersed spectra and the window 
filter was applied to avoid telemetry problems that could result 
from excessive counts in the zeroth order. A measure of pileup 
is the count rate per frame time (time to read out the chip) in 
a three-by-three pixel island. For this data set this quantity was 
always less than 0.02, and from Poisson statistics this implies a 
pileup rate of less than 1 per cent. The downside of using this 
mode is that spatial information is lost and the background is an 
order of magnitude higher. 

During spectral extraction regions, either side of the dis- 
persed spectra were used to create a background spectrum to 
be used in the analysis. In addition to instrumental background 
components this background was also used in the removal of 
components due to scattering. 

The observation lasted over four orbital cycles and was di- 
vided into 4 and 10 phase groups. For each group, HEG and 
MEG spectra with their responses were extracted and first-order 
spectra used in the present study. The spectral resolutions were 
0.023 A and 0.012 A for MEG and HEG, respectively (corre- 
sponding to 1 160 km/s and 580 km/s at the Si xiv line). 

The spectra were analysed with XSPEC (version 11) and ex- 
tracted into a format which could be used by programs written 
for IDL (version 6.2). We concentrated on the Sixiv 6.185 A 
Ly^-line as it is free of blends and the line profiles should be re- 
liable. However, profile fits for the S xvi Lya line at 4.73 A are 
also included in this work as well as the radial velocity study of 
the iron lines at 1.78 A and 1.86 A. 



3. XSTAR simulations 

XSTAR is a computer code that calculates radiation effects in 
a cloud surrounding an ionizing source ( Kallman 2()06j): in par- 
ticular, the complex emission line spectrum of different species 
is computed. The important input parameters are the chemical 
composition, density, and radius of the spherical cloud, and the 
X-ray/EUV spectrum of the central source. Of particular rele- 
vance is the radiation (ionizing luminosity) above the hydrogen 
ionization limit of 13.6 eV. The code was run in constant density 
mode, and the temperature was calculated to be the equilibrium 
value. The input parameters are listed in Table [T] and explained 
below. 



Table 1. The ionizing sources and gas parameters used in the 
XSTAR code. 



compact star 




WRstar 


gas envelope 


compPS 




BB 


N„=5x 10 23 cm" 2 


Tm,=0.71 keV 


T, 


,„ = 0.01 keV 


N=10" cm" 3 


L=2.46xl0 38 erg/s 


L: 


=5x 10 38 erg/s 


H=0.1,He=10,CNO 



Since the source of the wind is a WN(4-7) star (van 
Kerkwijk et al. 119961 we used helium-rich gas abundances 
with CNO equilibrium values: H=0.1, He=10, C=0.56, N=40, 
0=0.27; the other elements were set to solar, i.e. 1 (Crowther 
|20071 Vilhu 120071 ). The role of abundances is not significant, 
however, since we are not concerned with absolute values of 
fluxes, abundances or densities. The column (N# = 5xl0 23 
cm -2 ) and particle density (N=10 n cirT 3 ) were chosen result- 
ing in a cloud of size ~ 1 AU. The N^-value is compatible with 
the Thomson scattering origin of the orbital modulation of the 
hard X-ray light curves (Vilhu et al. 120071 1 . 

Both binary components were used as the ionizing radia- 
tors, thus forming a hybrid ionizer. To obtain initial estimates 
for the parameters of the X-ray ionizing source (i.e. the com- 
pact star), we fit RXTE/PCA spectra (xp5100, xp5200, xp5300) 
that were simultaneous to the Chandra observation in XSPEC 
with a Comptonizing model (compPS, Poutanen & Svensson 
1996) attenuated by photoelectric absorption (wabs). For a dis- 
tance of 9 kpc, the resulting unabsorbed bolometric luminosity 
is 2.46xl0 38 erg/s. The model includes photoelectric absorp- 
tion (8xl0 22 cm 2 ), a black body radiator (T bb = 0.71 keV), a 
Comptonized component (with kT e =5.8 keV and r=2.0) and a 
reflection component (with relative reflection = 4.5). The spec- 
trum is soft and dominated by a disc black body, typical for X- 
ray binaries in the high state. Note that we do not claim that this 
model is the best one to represent Cyg X-3 in the high state - the 
true absorption is probably more complex than is assumed here 
and the compPS model too simple. However, the model does 
give a viable numerical framework for our present purposes, in- 
cluding a correct value for the ionizing luminosity. 

For the WR star we used values similar to those of V444 
Cygni, which is a WN5 + 06 eclipsing binary: a mass of 10 Mq, 
a radius of 0.92 R , a luminosity of 10 5 1 L Q and T eff = 1 15 000 
K (Langer |19891 . We approximate the spectrum with a black 
body of T eff = 0.01 keV (1 16 000 K) and bolometric luminosity 
of 5xl0 38 erg/s; the ionizing luminosity is then 4.47xl0 38 erg/s. 
The incident and outgoing fluxes are illustrated in Figures Q]and 
|2] The outgoing flux was reddened by Ay = 20 and Ay/N# = 
5.3xl0~ 22 cm 2 , following Tables 21.6 and 21.7 in Cox (|2000b . 

Figures Q] and [2] show that even though the model outlined 
above adequately fits the observed Chandra HEG spectrum, it 
does not fully reproduce the infrared (IR) continuum and there is 
a soft X-ray excess visible relative to the continuum. Depending 
on whether Cyg X-3 is in quiescence or in outburst, the infrared 
magnitudes can vary by a magnitude (Hanson et al. 2000). FigQ] 
shows the infrared magnitudes given in Hanson et al. j2000| >: 
mH=13.1±0.1 and rns-=l 1.7+0.1, based on observations during 
quiescence by Fender et al. ( 1 19961 1. The soft X-ray excess is 
likely due to background issues and/or an inadequate representa- 
tion of the absorption in the RXTE/PCA spectra. However, these 
discrepancies do not affect the results presented here. 
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Fig. 1. The outgoing broadband spectrum. The dotted line traces 
the original spectrum and the solid line the reddened spectrum as 
computed with the hybrid ionizer (see the text and Table[TJ. The 
observed IR magnitudes (Hanson et al. 120001 short horizontal 
lines) and the mean Chandra HETG spectrum (OBSID 7268) are 
shown. The solid continuous line is the spectrum of the incident 
hybrid source. 




LOG (WAVELENGTH) ANGSTROEMS 



Fig. 2. A zoom of the simulation of FigQ] (shifted vertically 
by -0.4) overplotted with the mean Chandra HETG-spectrum 
(OBSID 7268). 



XSTAR requires as input a single ionizing source in the 
cloud center, although in our case we have two distinct ioniz- 
ers. Hence, several runs with different He-star luminosities and 
distances from the compact object were performed (using spher- 
ical shells) so as to interpolate the correct radiation field at any 
spatial point where the line emissivity was computed (see Fig. [3] 
for Si xiv). 

Additional emissivity grids were computed by varying the 
particle density. The degree of photoionization is governed by 
the density of the gas through the the ionization parameter, 
which is proportional to the ionizing luminosity and inversely 
proportional to distance squared and particle density, £=L/(nR 2 ). 
Ionization is sensitive to the density of the gas (which itself de- 
pends on the ionizing species in question): it can increase or de- 
crease in accordance. The final line emission was then computed 
in the wind model by fixing the mean density to 10 11 crrT 3 (see 




Fig. 3. Si xiv Lya-line emissivity per unit volume for constant 
density N= 10 1 1 cirT 2 (on orbital plane) based on the hybrid ion- 
izer (Table [T). The plot size is 32 R0, which corresponds to ten 
binary separations. The X-ray shadow of the WR-companion is 
seen on the right. The compact component is in the dark spot 
slightly to the left of figure center. In this grey scale white cor- 
responds to maximum and black to zero emissivity. At orbital 
phases 0.0 and 0.5 the system is viewed from the right and left, 
respectively. 



Fig. [4] for Si xiv). This mean density gives the correct electron 
scattering opacity along the line of sight and produces a partial 
eclipse in the light curve by a factor of two, as observed. A sim- 
ilar procedure was undertaken to estimate the absorption in the 
transmitted spectrum. 

The emission scales relative to the equilibrium temperature 
as T~° 6 , which explains the enhanced emission when the extra 
heater, the WR star, is farther away (see Fig. [3]). The line absorp- 
tion in the transmitted flux, in turn, has no temperature depen- 
dence. 



4. Particle simulations 

For modeling the velocity and density fields so as to obtain a 
reasonable input model for XSTAR, test particles were injected 
symmetrically from the Wolf-Rayet companion and their trajec- 
tories were followed in the binary potential. Neither acceleration 
nor interaction between the particles was assumed, which is a 
good assumption for the supersonic velocity of the wind. Typical 
values for the wind velocity and mass loss of WR stars are 2000 
km/s (observed wind velocities lie in the range 1000-3000 km/s 
(Cox 120001 )) and 3xl0~ 5 M /year, respectively (Langer [19891 . 
The masses of WN stars are between 10 - 80 Mq, while the 
range of WC stars is narrower (9-16 Mq, Crowther |2007l >. 
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Fig. 4. Si xiv Lya line emission in the photoionized wind model 
with q=5 and V=2000 km/s on the orbital plane. The plot size 
is 32 R0 (ten binary separations). The grey scale is linear with 
white corresponding to maximum and black to zero emission. 
At orbital phases 0.0 and 0.5 the system is viewed from the right 
and left, respectively. 



Table 2. Binary parameters used for wind particle simulations 
with wind velocities V = 1200 - 2200 km/s. 



M WR 


M c 


mass ratio q 


separation a 


V c 




(M©) 


CM©) 


(M WR /M C ) 


(R©) 


(km/s) 


(km/s) 


10 


2 


5 


3.2 


691.0 


138.2 


10 


5 


2 


3.4 


595.5 


297.7 



In our simulations we used wind velocities between 1200 
- 2200 km/s, a mass of 10 Mq for the WR star (with radius 
0.92R Q , Langer [1989j >, and 2 Mq and 5 M Q for the compact 
star (see Table [2]). We used 16-s time steps and the number of 
particles within the escape shell, 30xseparation, was kept con- 
stant at 10 6 . A new wind particle was always generated when 
a test particle moved outside the escape shell, returned back to 
the WR star, or passed closer than 0.15 Rq from the compact 
star (the assumed accretion radius). In all cases a stable structure 
was formed after about five binary revolutions. 

At large distances (larger than the binary orbit) the simulated 
outward directional velocity and particle density as a function 
of the distance, r, from the WR component, follow quite closely 
the semiempirical laws: 

V = const(l-R wr /r)^ ; N = const/[r 2 x(l-R wr /r/] ; y3=2. 



Table 3. Radial velocity parameters for the emission component 
of the Si xiv (6.185 A) and S xvi (4.733 A) Lya lines and mean 
^-values of P Cygni profile fits. 



line 


r 


K 


bluest 


MEG 


HEG 




km/s 


km/s 


phase 


* 2 (DOF) 


* 2 (DOF) 


Si xiv 
S xvi 


64+38 
75±33 


185±52 
100+36 


0.52±0.04 
0.44±0.08 


94(26) 
34(26) 


62(59) 
40(59) 



5. Line prof ile fitting results 

5.1. P Cygni model 

As already remarked in McCollough et al. (2000) and Liedahl et 
al. d2000T > the first HETG spectra of Cyg X -3 showed P Cygni- 
type line profiles (see also Sako et al.[T999). This is evident also 
from our phase-resolved Si xiv spectra. We modelled the line 
profile with two gaussians, one representing emission (gauss 1) 
and the other absorption (gauss2), as follows: 

F = exp(-gauss2)x(gaussl + continuum), where 

gauss = A0xexp-0.5x[(^-Al)/A2] 2 . 

The A-parameters are constants obtained in the fitting proce- 
dure and A is the wavelength. The continuum was kept constant 
over the line profile. In practice, A was replaced by the veloc- 
ity (8A/Axc) measured from the line center. This approach is an 
approximation of the P Cygni profile. The shape of the feature 
is a function of the distribution of wind velocities and is more 
complicated than the simple model above. 

The observed profiles were fit using the model described 
above with the IDL procedure mpcurvefit.pro (written by Craig 
Markwardt). Normal weighting was used (weights = 1/error 2 , 
where 'error' is the observation error of the flux inside a wave- 
length bin). The line profile was treated within + 4000 km/s 
from the line center, corresponding to 66 HEG and 33 MEG A- 
bins at the Si xiv line, respectively. In the final fit the absorption 
line central velocity was frozen to -900 km/s, to avoid artificial 
jumps in some parameters. However, this value is close to what 
is found when all seven parameters were left free. Mean intrinsic 
FWHM's of the absorption and emission profiles were 750+200 
km/s and 1850+200 km/s, respectively. Due to the requirement 
of high S/N, we were restricted to using only four phase bins. 

5.2. Radial velocities 
5.2.1. The Si xiv line 

The Si xiv Lya line profile fitting results are shown in Fig.|5]and 
Table [3] In Fig. 5 the dash-dot line shows the emission com- 
ponent (continuum + gauss 1) while the dotted line shows the 
absorption (continuum - gauss2). In Table [3] we also include 
the S xvi Lya-line. The x 1 values from the S xvi line fits are 
somewhat smaller than those of the Si xiv line - this is due to 
larger observational errors in the latter. Fig. [6] compares the 
observed Si xiv emission radial velocity curve with the model 
(q=5, V=2000 km/s, i = 60 deg). The rather conservative error 
bars were computed as errors in the sine curve fit to the mean 
MEG+HEG data (the solid curve). 
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-2 E , , , i , , , i , , , i , , , 3 

-4000 -2000 2000 4000 

VELOCITY (KM/S) 

Fig. 5. Si xiv 6.1850 A Lya line profiles at four orbital phases 
(HEG 1st order) overplotted with P Cygni-fits. The profiles are 
scaled with the maximum value and shifted vertically by a con- 
stant amount. The solid line gives the fit, dash-dot line the emis- 
sion component and the dotted line shows the absorption. The 
x-axis velocity is measured from the line centrum. 



600 



400 - 




-600 L , , , i , , i , , , , i , , i, 1 

0.0 0.5 1.0 1.5 2.0 

PHASE 

Fig. 6. Radial velocity curve of the Si xiv Lya emission compo- 
nent at four orbital phases overplotted with the model q=5 and 
V=2000 km/s for inclination 60 degrees (dashed line with er- 
ror bars). The data are represented by symbols with error bars: 
MEG=crosses, HEG=diamonds. 



The observed radial velocity semiamplitude of the Si xiv 
emission component was 185 ±52 km/s with a Gamma-velocity 
of 64+38 km/s and maximum blue-shift occurring at phase 0.5. 
The S xvi Lya line at 4.73 A behaves in a rather similar way but 
with a somewhat smaller radial velocity amplitude (see Table[3]). 

The maximum blueshift at phase 0.5 is due to the enhanced 
emissivity behind the X-ray source (Figs. [3] and @}. To explain 
the radial velocity amplitude observed, high wind velocity is re- 
quired. The effects of inclination, wind velocity and mass ratio 
are compared in Table |4j where ^ 2 -values between the observed 
and model radial velocity curves are computed for different in- 
put parameter values. Due to the conservative errors used, the 
^-values are relatively small. If the mean errors of the observa- 
tions were to be used, then the x 2 values would be larger by a 
factor of 5.2. 



Table 4. Effects of inclination, wind velocity and mass ratio on 
the observed and computed Si xiv radial velocity curves. 



inclination 




q 


2 

X 


degrees 


km/s 




DOF=28 


60 


2000 


5 


0.84 


45 


2000 


5 


3.7 


30 


2000 


5 


13.5 


60 


1800 


5 


19.9 


60 


2000 


2 


6.73 



Table 5. Observed (OBSID 7268 HEG) and model (q=5, 
V=2000 km/s) line luminosities (in units of 10 35 erg/s). 



ion 


wavelength (A) 


Observed Lum 


Model Lum 


Si xiv 


6.185 


1.00±0.08 


1.30 


S xvi 


4.733 


1.80±0.10 


1.60 


Fe xxv 


1.859 


7.50±0.60 


2.10 


Fe xxvi 


1.780 


6.05+0.50 


0.90 



5.2.2. The iron lines 

In the Fe-line complex at 6.6 - 6.9 KeV, the dominating H- and 
He-type lines are resolved in the HETG-spectrum (see Fig. [2). 
The lines are formed close to the compact star in a region of high 
ionizing flux and require special attention. In particular, our ion- 
izing wind models predict a much too weak Fe xxvi line, as can 
be seen in the large discrepancy between the observed line lumi- 
nosity and the model (see Table |5). In order to account for this 
discrepancy, an extra radiation source is required. The He-type 
Fexxv line is a composite of forbidden (1.869 A), intercombi- 
nation (1.859 A) and resonance (1.852 A) lines arising partially 
in the wind (with intensity ratios 0.35:0.45:0. 19) and partially in 
the region where the Fe xxvi line originates. 

To study the radial velocity curves we divided the observa- 
tion into four and ten phase bins. The line centroids and fluxes 
were determined by single Gaussian fitting in the HEG spec- 
trum. No absorption component was detected, since the spectral 
resolution is over three times worse than in the Si xiv line region. 
The wind model predicts for the Fe lines similar radial velocity 
curves to the cooler Si xiv and S xvi ions, but in fact the observed 
radial velocities don't correspond to the model. For Fe xxvi, the 
radial velocity analysis was performed separately for the four 
and the ten phase bins. In addition, we eliminated the two most 
blueshifted phase bins and undertook a third study without them 
(the eight bin case). The results are shown in Fig.[7]and Table|6] 

The maximum blueshift of the Fe xxvi line occurs around 
phase 0.25, when the compact star is approaching. Hence, it is 
tempting to associate the line-forming region with the compact 
component, presumably its accretion disc or corona. If the fit 
for Fe xxvi is forced to result in a maximum blueshift at exactly 
phase 0.25, the resulting K-amplitudes do not deviate signifi- 
cantly from those given in Table [6] 

Despite the apparent complexity in producing the Fexxvi 
line, we assume that it originates in the vicinity of the compact 
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Table 6. Parameters for the radial velocity curves of Fexxvi 
(1.78 A) and Fexxv (1.86 A) lines (treated as single gaussians 
and using different numbers of phase bins). 



line 


bins 


r 


K 


bluest 


redy 2 (DOF) 






km/s 


km/s 


phase 




Fe xxvi 


4 


28±95 


484±133 


0.27+0.04 


2.85(1) 


Fe xxvi 


10 


17±86 


454+108 


0.23±0.04 


0.96(7) 


Fe xxvi 


8 


102+92 


317+126 


0.24±0.04 


0.51(5) 


Fe xxv 


10 


324+56 


125±78 


0.10+0.09 


0.98(7) 




1 .0 



0.5 



0.0 I 

0.0 



0.5 



1.0 
PHASE 



2.0 



Fig. 8. The integrated area of the Si xiv-line absorption compo- 
nent is represented by crosses, the flux of the emission compo- 
nent by triangles. The diamonds joined by the dotted line gives 
the continuum at the Si xrv line. All parameters are scaled by 
their maximum values. 




< ; : 

< -2000 E , , , , , , , , , , , , , , , , , , , E 

0.0 0.5 1.0 1.5 2.0 

PHASE 



Fig. 7. Upper plot: Radial velocity curve of the Fe xxvi (H-like) 
Lyo'-line overplotted with the sine curve fit using 10 phase bins 
(crosses and the solid line). For comparison the same is shown 
when the observation is divided into 4 phase bins (sqaures and 
the dotted line). Lower plot: The 10 phase bin case for Fexxv 
(He-like). Lines were treated as single gaussians. Phases 0.0 and 
0.5 correspond to the cases when the WR star and compact star, 
respectively, are in front. 



object, primarily because the maximum blueshift occurs at phase 
0.25 when the compact object is approaching. 

For Fexxvi we compromise and adopt the mean K value 
in Table[6](K = 418 +123 km/s) to compute the mass function 
using the definition 

f = K 3 P/(2ttG) = M 3 xsin(i) 3 /M 2 of 

for a circular orbit. From this we get fwR = 1.51 ^1 ^O- ^ 
three cases for Fe xxvi lie, however, within the error ranges. To 
check this, we also analysed another high state observation of 
Cyg X-3 (OBSID 6601) that we divided into ten phase bins. The 
result for Fe xxvi was consistent with those of Table [6j K = 470 
+ 130 km/s when the maximum blueshift was frozen to phase 
0.25 (reduced^ 2 = 1.3). 

The He-like Fexxv line is more problematic. It is formed 
partially in the photoionized wind, like Si xiv, and in the vicinity 
of the compact object, like Fexxvi. However, extra complica- 
tions arise due to its triplet nature (resonance, intercombination 
and forbidden components) which may be selectively absorbed 
as a function of the orbital phase, and hence we do not treat it 
further here. 



5.3. Light curves of the emission and absorption components 

In Fig. [8] the integrated area of the Sixiv line absorption com- 
ponent (integrated gauss2) is represented by crosses (from the 
mean MEG+HEG data) scaled by its maximum value (525 
km/s). The continuum (scaled by 5.0 x 10 36 erg/s/keV) and the 
emission component (scaled by 1.5 x 10 35 erg/s) are represented 
by diamonds and triangles, respectively. 

Approximately one half of the Si xiv line emission originates 
from the far side of the X-ray source (compact star), accounting 
for the roughly similar orbital modulation of the line photons 
and the continuum. The wind particles themselves are relatively 
symmetrically distributed around the WR component (but asym- 
metrically around the X-ray component), thus causing the ob- 
served orbital modulation via electron scattering and photoelec- 
tric absorption. The Si xiv line absorption component is largest 
around phase 0.5 when the absorbing column is behind the com- 
pact star and not penetrating into the higher density wind region 
around the WR star. 

To study the iron lines (where no P Cygni analysis could be 
used) in the context of the Si xiv line, Fig. [9] was constructed 
using single gaussian fits with ten phase bins. As can be seen, 
the light curves for the lines are rather similar, suggesting the 
same absorbing agent in all three cases. This might be a com- 
bination of continuum opacity for line photons, such as electron 
scattering, plus broad line absorption probably in the underlying 
continuum originating from the X-ray source. 



6. Constraints on the system masses and 
inclination 

Hanson et al. (2000i) derive the mass function of the compact star 

f c = 0.027 +0.01 M = M c 3 xsin(i) 3 /(M c +Mi V/{ ) 2 , 

assuming that the infrared He i absorption line at 2.06 fim orig- 
inates in the base of the expanding wind (K = 109 +13 km/s) 
with expected maximum blueshift around phase 0.75, when the 
WR star is approaching). 
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Fig. 9. Scaled continuum fluxes (diamonds joined by dotted 
lines) and line fluxes (crosses, handled as single gaussians) as a 
function of the orbital phase for Si xiv (upper), Fe He-like (mid- 
dle) and Fe H-like (lower). Phases 0.0 and 0.5 correspond to the 
cases when the WR star and compact star are in front, respec- 
tively. 

This mass function can be combined with the mass function 
we derived from the Fe xxvi line above: 

f WR = 1.51^| M = M H , s 3 xsin(i) 3 /(M c +M lVR ) 2 . 

For the mass ratio we obtain Mwr/Mc = 3.8^'T. Constraints 
for the compact star mass and inclination are approximately 

M c = 2.8+!;5X(M WR /10) M Q and 
i = (38 +12) x(M ws /10) a33 degrees. 

More mass increases orbital motion; hence a massive system 
requires smaller inclination to account for the radial velocity am- 
plitude. As an example, if Mwr = 60 Mq then a possible solu- 
tion (within the error limits) can be Mc = 25 Mq with a small 
inclination of 15 degrees. This is illustrated in Fig. [TOl where 
the permitted (Mc, inclination) ranges are shown for three WR 
masses. 

Inclinations as small as 30 degrees are acceptable for charac- 
terizing the amplitude of the Si xiv line. However, the inclination 
can not be constrained in this way since smaller inclination can 
be compensated for by larger wind velocity (see Table @). For 
large wind velocities the shape of the continuum light curves is 
not very sensitive on inclination either. The effect of mass ratio 
on light curves is almost negligible at high wind velocities. 

A small inclination is supported by population synthesis re- 
sults on core He-burning wind-fed He+BH binaries by Lommen 
et al. (120051 1. They estimate that there is one such system in the 
Galaxy with an orbital period close to that of Cyg X-3 and sat- 
isfying the disc-formation criterion ( 1 1 Mq + 5 Mq, compare 
with our Fig. ITOt . The existence of a disc in Cyg X-3 is sup- 
ported by the observed superluminal jets. Lommen et al. ( 120051 1 
consider unlikely the existence of similar WR+NS binaries with 
discs, as already shown by Ergma and Yungelson (1998), due 
to the large wind velocity and small angular momentum of the 
WR+NS wind that hinder disc formation. On the other hand, 
Marti et al. (120011 1 detected a bipolar jet flow at a large angle 
(73 degrees) to the line of sight. If the flow is parallel to the 



Fig. 10. Permitted (compact star mass, inclination)-ranges for 
three WR-companion masses marked (in Mq) using mass func- 
tions from the Fe xxvi (this study) and He i (Hanson et al. 120001 ) 
lines. 

orbital plane axis, then this favours large inclination. However, 
this need not be the case for a wind-disc scenario as discussed 
by Marti et al. (2001) who assume an inclination of 30 degrees 
(see their Fig. 7). If this is indeed the case, then small inclination 
is favoured. 

Gladstone et al. ( 120071 1 consider a black hole signature to be 
a source penetrating into the ultrasoft domain in their colour- 
colour diagram (where Cyg X-3 never entered). However, if the 
mass of the black hole is small, then its disc is hot, keeping it out 
of the ultrasoft domain populated by more massive black holes. 
On the other hand, based on spectral evolution vs. the Eddington 
luminosity, Hjalmarsdotter et al. ( 12008bt point to a massive (30 
Mq) black hole, and consequently a very small inclination, with 
a massive WR companion. 

Support for low inclination in our photoionized wind model 
comes from the infrared. We computed the emissivity and radial 
velocities of the IR Hen 2.09 fim emission line (using atomic 
data provided by T. Kallman, private communication). It turns 
out that most of the line emission comes from the X-ray shadow 
(see Fig. 01 where the WR star is the only ionizer. The radial ve- 
locity K-amplitude is around 500 km/s with maximum blueshift 
at phase 0.0, i.e. 180 degrees from that of the Si xiv emission, as 
observed by Hanson et al. (12000b . However, in our model with a 
wind velocity of 2000 km/s, q=5 and inclination 60 degrees, the 
computed K value is much larger than that observed (1020 km/s 
with maximum blueshift at phase 0.0). If the inclination is 30 
degrees, the resulting amplitude of 580 km/s would be closer to 
that observed (see Fig. ITTt and the permitted ranges of the WR 
and compact star masses are 10 - 40 Mq and 2.8 - 8 Mq, re- 
spectively. For an inclination of 60 degrees the same ranges are 
5 -8 M and 1.0- 3.2 Mq. 

7. Discussion 

Our mass estimates (see Fig. ITOb corroborate the Hanson et al. 
(2000) conclusion that if the compact component is a neutron 
star with inclination greater than 60 degrees, then the WR mass 
must be in the range 5-11 Mq. Stark and Saia (2003]) gave an 
upper limit of 3.6 Mq for the compact star mass using the He- 
like Fe xxv line with a similar radial velocity curve to ours (see 
Fig. 7) - however, we prefer to use the H-like Fe xxvi line here 



8 



O. Vilhu et al.: Orbital modulation of X-ray emission lines in Cygnus X-3. 




Fig. 11. Model (q=5, V=2000 km/s) radial velocity curve of the 
infrared Hen 2.09 fim line for inclination 30 degrees. The line 
is formed in the X-ray shadow cone behind the WR component 
where the wind is photoionized by the EUV radiation from the 
WR star. The fits for the observed Si xiv and Fe xxvi (H like) 
lines are shown for comparison. These lines are formed in the 
X-ray photoionized wind and at the compact star, respectively. 



and defer a more in-depth discussion of the He-like Fe xxv line 
to a separate paper. They also observed radial velocity curves 
of the Si xiv and S xvi lines very similar to ours. The difference 
is in their larger Gamma (recession) velocities which are due to 
ignoring the P Cygni absorption and treating the lines as pure 
Gaussians. 

The formation mechanism of the Fe xxvi line at 1 .78 A re- 
mains unknown, although it was proposed above to arise in the 
vicinity of the compact object. If the line is formed in a hot col- 
lisionally dominated plasma with temperature T = 10 KeV, the 
required emission measure EM = n 2 V is around 10 60 cm -3 (for 
a distance of 9 kpc) using the MEKAL model in XSPEC. The 
emitting plasma can be e.g. a spherical corona with radius 10 10 
cm (roughly the assumed accretion radius) and electron density 
5xl0 14 cm 4 . The corona itself can also be a composite of many 
impulsive flares. Another alternative is that the hard coronal flux 
is photoionizing the disc providing additional Fe emission. 

In our particle simulations we assumed Mwr = 10 M0 and 
one may ask how a larger WR mass would influence the model 
results. A larger mass increases the separation and, hence, low- 
ers the wind density in line forming regions. Keeping all other 
parameters fixed, the major effect of this would be the decreasing 
of model line luminosities in Table [5] Unfortunately, the masses 
of WR stars have a large spread even for the same spectral type 
(see Fig. 4 of Crowther l20"0"7l 

The question of the origin of line absorption (Section 5.3) is 
interesting, even though not fully solved here. The similarity of 
the Si xiv, Fe xxv and Fe xxvi light curves (Fig. |9]l is indicative 
of the same absorption agent for all. The P Cygni treatment of 
Si xiv (Fig. [8} leads us to speculate that line photons are ab- 
sorbed in the same fashion as the continuum. In addition, there 
can be line absorption operating in the continuum (below the line 
photons) that is blueshifted and varies with the orbital phase. A 
combination of these two absorptions might then cause the broad 
minima seen in Fig. [9] We leave a deeper discussion on line ab- 
sorption for a separate paper. 



8. Conclusions 

We have shown that the Si xiv (6.185 A) and S xvi (4.73 A) Lya 
emission lines of Cyg X-3, using Chandra HETG observations, 
have P Cygni-type orbitally modulated line profiles. The emis- 
sion components arise in the wind of the WR-companion, pho- 
toionized and excited by a hybrid source, which we take to be 
composed of the compact X-ray source and WR He-star. 

The radial velocity curves of the emission components have 
maximum blue-shifts at phase 0.5, when the compact star is in 
front, as predicted by the wind model. The K-amplitudes depend 
mostly on the wind velocity and constraints on the mass ratio and 
orbital inclination are weak. The absorption component veloci- 
ties are constant around -900 km/s. Both the absorption and the 
emission components have intensity maxima (integrated gaus- 
sians) around phase 0.5, similar to the continuum flux. 

The Fe xxvi Lya line at 1 .78 A appears to arise in the vicinity 
of the compact star, as maximum blue shift occurs at phase 0.25 
which is when the compact star is approaching. Combining the 
mass function derived from this line with the mass function from 
the IR He 1 absorption line (Hanson et al. |2000| > we were able to 
constrain the system masses and inclination. Our photoioniza- 
tion simulations of the IR He 11 2.09 fim emission line (formed 
in the X-ray shadow) lend support to a small orbital inclination 
(30 degrees) for which the range of masses of the compact star 
is 2.8 - 8.0 Mq For a larger inclination of 60 degrees this range 
would be 1.0- 3.2 M©. 

In the light of our results, we either have a neutron star and 
low-mass Wolf-Rayet star in a system of large orbital inclination 
or a black hole whose mass depends on the Wolf-Rayet mass and 
the system orbital inclination. 
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